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ABSTRACT 

We present the first Spitzer Infrared Spectrograpli (IRS) observations of the 
[O IV] 25.89 /im emission line detected from the ultraluminous X-ray source 
(ULX) in Holmberg II. This line is a well established signature of high ioniza- 
tion, usually associated with AGN. Its detection suggests that the ULX has a 
strong impact on the surrounding gas. A Spitzer high resolution spectral map 
shows that the [O IV] is coincident with the X-ray position of the ULX. The ra- 
tios of the [O IV] to lower ionization lines are similar to those observed in AGN, 
suggesting that a strong UV and X-ray source is responsible for the photoion- 
ization. The best XMM-Newton data is used to model the X-ray band which is 
then extrapolated into the UV. We perform infrared and ultraviolet photometry, 
and use previously published optical and radio data to construct the full SED 
for the ULX and its companion. The preferred model to describe the SED in- 
cludes an accretion disk which dominates the soft X-rays but contributes little 
at UV and optical wavelengths. The optical counterpart is consistent with a B 
supergiant as previously suggested in other studies. The bolometric luminosity 
of the ULX suggests the presence of an intermediate-mass black hole with mass 
>85 Mq for sub-Eddington accretion or, alternatively, a stellar-mass black hole 
that is accreting at super-Eddington rates. In a follow-up second paper we per- 
form detailed photoionization modeling of the infrared lines in order to constrain 
the bolometric luminosity of the ULX. 



-2- 



Subject headings: black hole physics — galaxies: individual (Holmberg II) — 
infrared: ISM — X-rays: binaries 



INTRODUCTION 



The majority of our knowledge of ultraluminous X-ray sources (U LXs) stems from ove r 



200 observations of these objects by both Chandra and XMM-Newton (iLiu &: Mirabelll2005l ). 



The elevated X-ray luminosities of ULXs have raised intriguing questions about their nature 



19991): however substan t iating evidence for such objects is lacking (IStobbart et al 



Colbert & Musi 


lotzkv 


(Stobbart et al. 


20061; 



Goncalves fc Soriall2006l : iBerghea et al.l 120081). Attenti on has switched to less exotic sources 
for the emission, such as beaming (King et al.M200ll ) or super-Eddington accretion from 



a stellar mass black hole (IBegelmanl |2002| ). both of which explain ULX X-ray properties 
without an IMBH. The idea is not far fetched, since ex amples of super-Edding ton sources 
have been found in our Galaxy, such as GRS 1915-1-105 (iFender fc Bellonil 120041 ). V4641 Sgr 
( iRevnivtsev et al.ll2002l ) and possibly SS 433. Indeed SS 433 could be an example of both 
beaming and super-Eddin gton accretion, the combination of which cou ld easily explain even 
the most luminous ULXs (iBegelman et al.ll2006l : iPoutanen et al.l 120071 ). 



Very little work has been done on ULXs at non X-ray wavelengths, especially in com- 
parison with AGN and Galactic stellar-mass black holes. This is in part because, unlike 
standard AGN, the emission from the accretion process itself is harder to detect at these 
wavelengths. Indeed, the X-ray to optic al luminosity ratios in ULXs are at least ten times 
higher compared to those seen in AGN (jWinter et al.l 120071 ). To make matters worse, there 
are no ULXs in our own Galaxy (though SS 433 is regarded as a ULX by some authors), 
and the stellar ULX companions are very difficult to identify in other galaxies. 

What little multi-wavelength work done on ULXs thus far has been executed primarily 
with high resolution optical spectroscopy. In a few cases, ULXs have been found to have 
an observable effect on their local environment - developing ionized bubble nebulae of up 



to 500 pc diameter (jPakuU fc Mirionill2002l ). These are much larger than regular supernova 
remnants and energy estimates show that the bubbles are consistent with being blown either 
by relativistic jets from a ULX, or abnormally powerful supernovae. In some cases high 



excit ation lines have been detected, such as He II A4686 (LP. = 54.4 eV; e.g.lAbolmasov et al. 



20071 ). Only very hot Wolf-Ra yet stars can produce t his line (ISchaererlll996l ). but the ULX 
companions are not such stars (lAbolmasov et al.ll2007l ). Therefore, the detection of the He II 
line must be due to the interaction of the ULX with the surrounding medium. Interestingly, 
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SS 433 is also surrounded by a 100 pc radio nebula, W 50 (e.g. lFabrikall2004l ) 



The importance of ionization nebulae for investigating X-ray sources has long been 
reali zed, beginning with the fir st detection of nebulae around super soft sources such as LMC 
X-1 JPakuU fc AngebaultlEoseh and CAL 83 jRemillard et al.lll995h in the Large Magellanic 
Cloud (LMC). The high-excitation lines observed in nebulae, such as the He II recombination 
line, can act as giant bolometers for the ionizing UV and X-ray emission, thus allowing 
the intrinsic hard spectral energy distribution (SED) to be inferred. Long-term averaged 
luminosities can be estimated and the UV emission constrained. This provides tremendous 
value when studying the nature of ambiguous ULXs, because the luminosities estimated from 
ionization are independent of X-ray spectral modeling estimates. In fact, such luminosities 
are more accurate, because the X-ray-based estimates are in general lower bounds and depend 
on the viewing angle. More importantly, the ionization-estimated luminosities are time- 
averaged estimates, which is critical for variable sources such as ULXs. This is especially true 
for ULXs with very few X-ray observations, since some of these sources can be transient. On 
the other hand, photoionization modeling is dependent on many parameters which are often 
not known very precisely, such as the gas density, geometry, metallicity etc. Contamination 
from other ionizing sources is also a concern especially in crowded regions. However, when 
used together. X-ray analysis and photoionization modeling can provide good constraints 
on the ULX luminosity. For example, a simple comparison between the two estimates can 
constrain the beaming factor. 

One of the most interesting examples of an ionized nebula associated with a ULX 



is located in the dwar f galaxy Holmberg II (jPakuU fc Mirionil 120021 : iLehmann et al.l 12005 



Abolmasov et al.l 120071 ). He I I is detected in a s mall nebular region around the ULX, the 
"Heel" of the "Foot Nebula" JXaaret et al.lbooi KWZ hereafter). An OB star is the likely 
counterpart (KWZ). The morphology of He II emission and other optical lines is consistent 
with the ph otoionization of the nebula by the ULX. KWZ confirmed this with CLOUDY 
simulations. iMiller et al.l (120051 ) detected extended radio emission from a ~50 pc diameter 
region coincident with the position of the ULX. The authors conclude that the flux and mor- 
phology of the radio emission is inconsistent with emission from a SNR or an H II region, 
and is probably associated with the ULX activity. 

In this paper and a subsequent companion paper we use IR Spitzer spectroscopic ob- 
servations for the first time to constrain the underlying hard SED of a ULX. The IR ob- 
servations, when combined with previously published optical and X-ray data, suggest that 
the powerful UV and soft X-ray emission from the Holmberg II ULX photoionizes the sur- 
rounding medium. In Section 2 we present the Spitzer IRS data analysis and XMM-Newton 
X-ray spectral fitting results. We perform additional photometry measurements with the 
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Infrared Array Camera (IRAC), the Multiband Imaging Photometer for Spitzer (MIPS), 
the European Photon Imaging Camera (EPIC) and Optical Monitor (OM) instruments on 
XMM-Newton, and the Galaxy Evolution Explorer [GALEX). In Section 3 we analyze the 
possible contribution of self-irradiated accretion disk to the UV and optical fluxes. In Sec- 
tion 4 we use our data and other published data to construct an SED of the Holmberg II 
ULX. This is also used in the companion paper to further constrain the X-ray spectrum. In 
Sections 5 and 6 we discuss the significance of the [O IV] line detection, compare the IR line 
ratios with those of AGN and present the final observational results. A follow-up paper will 
present detailed photoionization analysis of the Spitzer data. 



OBSERVATIONS AND DATA ANALYSIS 



We begin this section with a description of the Spitzer IRS (IHouck et al.ll2004l ) obser- 
vation of the Holmberg II ULX. We next present multi-wavelength photometry and X-ray 
spectroscopy used in the subsequent section to construct the spectral energy distribution 
(SED) for the ULX. 



2.1. S'pifeer IRS Observations 



The data presented here are archival Spitzer data, taken on 2004 October 4, as part 
of the Legacy Program Spitzer Infrared Nearby Galaxies Survey (SINGS Program 193, see 
Kennicutt et al.ll2003l ). The IRS spectral map is described by three parallel pointings and 
five perpendicular at half-slit steps, for both Short High (SH, in the range 9.9-19.6 fim) and 
Long High (LH, in the range 18.7 — 37.2 /im) modules. The total exposure time was 240 s 
for both SH and LH observations. The data were preprocessed by the IRS pipeline (version 
11) at the Spitzer Science Center prior to download. Spectral maps were constructed using 
"BCD-level" proces sed data in conjunction with CUBISM (CUbe Builder for IRS Spectra 
Maps) version 1.50 (jSmith et al. 2007; Kennicutt et al. 2003 ). The absolute calibration error 
on fluxes obtained with CUBISM is ~25% JSmith et al.lbooi ). 



We adopt a distance of 3.05 Mpc to Holmberg II (KWZ, 1" is 15 pc). The LH maps 
obtained with CUBISM have intrinsic sizes of 44.6" by 28.8" (670 pc by 430 pc), while the 
SH maps are 27" by 13.6" (405 pc by 205 pc). However, in order to avoid aperture effects, 
the spectra presented here were extracted from like-sized apertures. Thus the SH and LH 



Available at: 



littp://ssc. spitzer. caltech.edu 



/archanaly/contributed /cubism / 
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spectral data result from an equally sized aperture of 23" by 13.6" (340 pc by 205 pc), which 
is only slightly smaller than the full SH aperture. We also extract spectra from the smallest 
aperture permitted by CUBISM: an 8.9" by 8.9" aperture centered on the ULX (i.e. 2x2 
LH pixels, or 130 pc by 130 pc). We will refer to these apertures as the standard aperture 
and the small aperture, respectively. 

The high-resolution spectra (10 /xm— 37 ^m) for the two apertures are plotted in Fig- 
ure dJ Since the spectra are extracted from matched apertures, the spectral lines in both 
the SH and LH modules are properly calibrated and the continuum is properly matched. In 
Figure [2] we show in detail the 4 detected lines, with the theoretic al line centers. We used 



SMART (Spectroscopic Modeling Analysis and Reduction Tool, see lHigdon et al.ll2004l ) ver- 
sion 6.2.6 to measure line fluxes. Only 4 lines were detected, the fluxes and upper limits 
for non-detections are presented in Table [TJ Here we show the results for both the standard 
aperture and the small aperture. We note that the [O IV] line is detected at a sig nal-to-noise 



ratio > 10. The spectral lines for the Holmberg II galaxy are also presented in iDale et al. 



( I2OO8I ). The focus of that paper was not the ULX in Holmberg II, and as a result, the 
extraction region for the relevant spectral lines is larger than the apertures presented here. 
However, wi th this in mind, we flnd that our line fluxes agree reasonably well with those 



presented in iDale et al.l (120081 ) . 



The SH and LH maps for the 4 detected lines are shown in Figure [31 A false-color 
Hubble Space Telescope (HST) image shows the Foot nebula in Ha emission, together with 
the standard and small extraction apertures. Because of the poor spatial resolution of the 
Spitzer spectral maps, the emission for most lines is difficult to disentangle from the nearby 
star-forming regions in the Foot nebula. This is especially true of the low ionization lines 
([Ne III], [S HI], and [Si II]). The morphologies of these maps indicate that the emission is 
likely coming from both the ULX and the surrounding star forming regions. However, while 
the low ionization lines [Si II] and [S HI] are very extended, the [O IV] emission is visibly 
compact. The peak of the of the [01 V] emission is coincident with the Chandra position for 
the ULX, with the majority of the emission concentrated in the 2x2 pixel region, overlapping 
with the X-ray point source. There is a slight extension to the west and south that is visible 
on the spectral map and which may imply that the [OIV] emission is extended. These faint 
extensions correspond very well with two regions of young star formation visible in the HST 
image in Figure [31 

KWZ found that the He II line, which has an ionization potential just slightly lower 
than [OIV] (54.4 eV, compared to 54.9 eV), is also concentrated in a small region of 50 pc 
diameter around the identified optical counterpart of the ULX. Because the spatial resolution 
is limited, a detailed comparison of our maps with the optical maps presented by KWZ cannot 
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be made. However, we find that the morphology of our [O IV] map appears similar to He II, 
being brighter in the east and fainter and more extended in the west. 



2.2. X-ray Observations and Spectral Fitting 



We use X-ray spectroscopy to model the X-ray emission from the ULX. We also use it to 
estimate the UV emission by extrapolating the model to lower energies. X-ray observations 
of the Holmberg II galaxy are available from ROSAT, ASCA, XMM-Newton and Chandr a 
jZezas et al.1ll999l : iMivaji et allbnOll : IXerp et allbnod : IXaaret et ahlbond : lOoad et allbnOBl ) . 



The X-ray luminosity of the ULX is in excess of 10^° erg s~^ in all observations, with the 
exception of a peculi ar low and soft state i n September 2002, when the luminosity dropped 
to ~5xl0^^ erg s~^ ( jPewangan et al.l 12004 ). 



To model the X-ray spectrum we used the best available data from the longest (~100 ks) 
XMM-Newton EPIC PN observation in th e archive (Ob s ID 0200470101, taken on 2003 
December 5). The data were processed by IWinter et al.l (120061 ). We used XSPEC version 
12.4 to model the spectrum between 0.3 and 8.0 keV. The X-ray spectrum is the most 
important component of the total SED since it is the only observed emission that directly 
emanates from the ULX, and is not contaminated by the surrounding star-forming region. 
For this reason we have taken many precautions in modeling the X-ray spectrum. 

The (far) UV emission from the ULX is the least known, but the most important for 
ionization. The only available method to estimate the shape of the UV is to extrapolate 
from X-ray models. Therefore, it is very important to find the right X-ray model to fit the 
data, especially in the soft X-ray band. Three physically different models were used. Each 
model contains different information about the accretion process. The first model is very 
basic and is widely used to model the accretion process for stellar mass black holes, ULXs 



and other accreting objects (IMiller et al.ll2004l : iMcClintock fc Remillardll2006l : IWinter et al. 



2006 



Feng &: Kaaretl l2006l ). The second is a modification of the first model that is more 
physically realistic when extrapolated to lower energies. Finally, the third model is more 
flexible than the other two models in the UV and allows for a tighter fit in the UV while 
staying within the X-ray error bars. The intrinsic (unabsorbed) X-ray luminosity is >10^° 
erg s^^ for all models. The results from the X-ray spectral fits are listed in Table [21 and a 
summary of the models is as follows: 

1. Power-law Plus Multicolor Disk (PLMCD) Model: This is a tw o-component 
model with a multi-color accretion disk component (MCD, Mitsuda. et al.lE984 ) and a power- 
law for the hard tail. Following KWZ, the absorption column density was fixed at the Galac- 
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tic value of 3.42x10 cm~ , and an addit ional local absorption (N^) model was a dded. We 
assume a local metallicity of 0.1 Zq (see iKaaret et al.ll2004l : iPilyugin et al.ll2004l ). This is 
a very simplistic model commonly used to fit X-ray data for ULXs. However, the model 
falls apart when extrapolated to lower energies, where the power-law dominates. The ex- 
trapolation is actually physically u nrealistic, since it i s imp ossible for the disk component to 
Comptonize lower energy photons. iDone fc Kubotal (120061 ) pointed out that the power-law 
component in such cases needs to have a break at low energies, in agreement with the theo- 
retical model where the photons from the disk are up-scattered in the hot corona. Moreover, 
if the power-law is produced by another mechanism (e.g. a jet), it would not extend down 
to the soft UV and optical spectrum because the fluxes here would be over predicted (see 
below. Section 4). Nevertheless, the model continues to be used extensively in X-ray spectral 
fitting of ULXs. In conclusion this model does fit the X-ray data well (Ax^/dof = 1.14), but 
it does not make physical sense when extrapolated to lower energies. It is shown here as a 
base, from which alterations can be made, resulting in more physically meaningful models. 
The following model is one such modification to the PLMCD model, and is more physically 
realistic when extrapolated into the UV. 

2. Modified PLMCD Model: In order to isolate the intrinsic disk emission from the 
Comptonized disk emission, we fitted the MCD and power-law separately for two different 
regions of the X-ray spectrum. We first fitted the hard tail with a power-law. The spec- 
trum departs from this power-law around ~1 keV. Then we fitted an MCD model to the 
0.3—1.0 keV band. The disk temperature (0.38 keV) is slightly higher than other published 
fits (KWZ obtained ~0.2 keV), because our fitting methods are slightly different. Lower 
temperature accretion disks are generally more effective ionizing sources. Therefore this 
higher temperature model may be viewed as a lower bound to the SED in the FUV and 
soft-Xray bands. We note that the Modified PLMCD differs from the PLMCD model in 
that the UV band is predicted from the disk component rather than the extrapolated PL 
component. This difference yields a much more physically realistic interpretation for the 
emission. 

3. Broken Power-Law (BP L) Model: Finally w e also fitted the X-ray spectrum to 
a simple broken power-law model. iKaaret et al.l (120031 ) fit a similar model to the ULX in 
NGC 5408. Such a m odel could be explained if th e emission from the companion star were 
Comptonized by a jet (jCeorganopoulos et al.ll2002l ). This model is much more flexible in the 
UV than either the PLMCD or the Modified PLMCD models, and can predict the slope of 
the SED in the UV while simultaneously fitting the hard power-law component in the X-rays 
within the error bars of the X-ray measurements. Surprisingly, this simple model provides 
a very good fit to the X-ray data, better than the PLMCD model. The break energy is at 
1.0 keV, the same energy at which we split the two components in the Modified PLMCD 
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model (2 above). Finally, we note that, even if this model has a power-law less steep than the 
PLMCD model at lower energies, it is still overpredicting the soft UV and optical emission 
according to our photometry results (Section 4, below ) . A s imilar overprediction at lower 
energies was found for this model by iDone &: Kubotal (120061 ) for the Galactic binary XTE 
J1550-564 when at high accretion rates (in the very high state). 



2.3. SED of the ULX: Multi-wavelength Photometry 

In order to constrain ULX multi- wavelength SED, we performed photometry in the IR 
and UV using Spitzer, XMM-Newton and GALEX data. We also incorporated previously 
published optical and radio data from HST and the Very Large Array (VLA). A short 
description of the data and the measurements are shown in Table [31 All our photometric 
measurements were extracted from an aperture of 50 pc diameter (3.3"). We used the 
optical counterpart identified by KWZ to center the extraction aperture (see Fig. [3]) . Aspect 
correction for the IR and UV images was performed using bright sources in field, such as 
the foreground star seen in Figure [31 The measurements presented in Table [31 represent 
the integrated fluxes over this region, including the nebular emission and the contamination 
from other nearby sources. They should therefore be considered upper bounds to the intrinsic 
ULX contribution. 

IR Photometry: IRAC and MIPS images (see Fig. 4) were taken as a part of SINGS 
Program 195. We extracted 3.6 //m and 8.0 /im fluxes using 50 pc diameter apertures centered 
on the ULX position. For the MIPS 24 /dm photometry, we used a 100 pc extraction aperture 
and then divided the flux by four to obtain an approximate flux for an equivalent 50 pc 
aperture. We followed the standard extraction procedure and conversion factors provided 
in the IRAC Data Handbook^. The pipeline mosaics for IRAC and MIPS are accurate to 
within 20%. 

UV Photometry: XMM-Newton Optical Monitor (OM) images of Holmberg II (Fig. 4) 
are available in the UVWl band, which is centered at 2910A(F WHM 620A) We use the 



mosaicked image created for the Optical Monitor Catalog (OMCat. lKuntz et al.ll2008l ). which 
has an improved coordinate correction (better than 0.5") over the automated pipeline. To 
measure the UVWl count rate we used the XIMAGE software included in the HEAsoft 
package. We converted the count rate into flux using the conversion factors listed in the 



^http://ssc. spitzer. caltech.edu/irac/dh/ 
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XMM-Newton Science Analysis System User's Guid^. We obtained a UVWl magnitude of 
18.42 ± 0.06, which corresponds to a luminosity of 8.3 xlO'^'' erg s~^ and a flux density of 
33.1 ± 0.1 /iJy. The OM image in Figure 4 shows that the Foot Nebula has a strong UV 
emission which peaks on the OB association of stars seen in the HST image of Figure [3] west 
of the ULX. The total UVWl flux for this UV bright region (including the ULX and its 
optical counterpart) is equivalent to 29 05V stars. Figure [3] shows that these bright stars 
have apparently blown a cavity in the NW part of the nebula, thus producing the "Sole" 
and the "Toes" of the Foot nebula. We note that these cavities are correlated with the 
faint extended [O IV] emission seen in the CUBISM spectral maps (Fig. [3]). The strongest 
emsission of Ha is however in the Heel that hosts the ULX. 

We also used archi ved GALEX observ ations, taken as a part of the GALEX Nearby 
Galaxies Survey (NGS, iBianchi et al.ll2006l ). Photometric fluxes were extracted from the 
FUV (1350-2013 A) and NUV (1750-2013 A) bands with the XIMAGE too\. The GALEX 
images have a lower resolution compared to OM, but they show similar features. 

Published Radio Data: Finally, we searched the liter ature for pub l ished photomet 



ric data for the Holmberg II ULX. Using radio VLA data, iMiller et al. 



J2005h measured 



1.174 ± 0.0085 mJy at 1.4 GHz and 0.677 ± 0.207 mJy at 4.86 GHz. Miller et all (120051 ) 
found extended emission covering the Ghandra position of the ULX: 3.7" by 2.7" at 1.4 GHz 
and 3.4" by 1.9" at 4.86 GHz. These are comparable to the 50 pc extraction region used in 
the IR and UV photometry. 

Published Optical Photometry: In addition, KWZ measured a V-band magnitude 
of 22.04 ± 0.08 for the stellar optical counterpart from observations with the Advanced 
Camera for Surveys (ACS, F550M filter, see Fig. [3]) on the HST. The optical flux is likely 
less contaminated by the nebular emission or emission from other sources compared to our 
measurements, because of the much smaller aperture. 

Extinction Corrections: The Milky Way reddening in the direction of Holmberg II is 
E(B— V) = 0.032 (from the NASA Extragalactic Database). KWZ estimated an additional 
local reddening of E(B— V) = 0.07 ± 0.01 based on the X-ray absorption. Thus the total 
extinction for Hy = 3.1 is Ay = 0.32. We show the projected extinction and the extinction- 
corrected values for the U V and optical data in the last two columns of table Table [31 These 
were calculate d using the ICardelli et al.l (119891 ) extinction curves. For the infrared data, we 
estimate from iDraind (119891 ) that the extinction is A^ < 0.02, and therefore negligible. We 
note that the photometric fluxes presented here are likely contaminated by the surrounding 
nebular emission as is evidenced by the HST images in Figure [31 Moreover, KWZ caution 



•^http: / /xmm2.esac.esa.int/external/xmm_user .support/documentation/ 
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that the local extinction could in fact occur very close to the compact object and therefore not 
affect our measurements at all. We will therefore take the (extinction-corrected) photometric 
measurements as upper bounds. 



3. SELF-IRRADIATION OF THE ACCRETION DISK 

While the MCD accretion disk model is a good approximation of the UV continuum to 
first order, it does not account for self-irradiation of X-ray photons caused by thickening of 
the accretion disk at the outer edges and relativistic effects on the X-rays emitted from the 
inner regions of the disk. In Galactic X-ray binaries, this reflected radiation has been shown 



to increase the optical luminosity of the accreti on disk by as much as 18% (ISanbuichi et al. 



19931 : IVrtilek et al.lll990l : iRevnivtsev et al.ll2002l ). Thus if irradiation were contaminating the 



UV emission in the Holmberg II ULX, such a model must be included in the final fit. 

We therefore checked to see if X-ray irradiation would affect the ionizing UV emission 
from the Holmberg II ULX , or co ntribute significantly to the flux of the optical counterpart. 



Following ISanbuichi et al.l (Il993l ). we modeled the expected SED of an irradiated disk for 
black holes of mass 10, 100 and 1000 Mq and with varying accretion disk radii. Here the 
accretion disk radius is constrained by the size of the Roche lobe, and is therefore dependent 
on the spectral type of the companion star. We find that for stellar mass black holes, 
irradiation affects the accretion disk model significantly for most companion spectral types. 
However, for intermediate black hole masses, irradiation is only important to systems with 
late-type supergiant companions. 

In Figure 5 two examples of irradiated spectra are plotted. These represent typical 
spectra for stellar-mass BHs (10 Mq) and IMBHs (1000 Mq). In the latter case, an accretion 
disk typical for a system with an M-supergiant was used to illustrate the effect of irradiation. 
Only if the IMBH has such companions is the self- irradiation effect significant. However, 
we note that an M-supergiant companion is not consistent with the observations for the 
ULX in Holmberg II (see Section 4 below). As can be seen from this plot, in both cases 
the irradiation dominates the optical and near-IR spectrum much more than the UV band. 
Thus we conclude that while irradiation may be important to model in the optical spectrum 
of a system with a stellar mass black hole, it is likely not dominant for an IMBH system. 
Finally we conclude that the irradiation does not affect the detected IR lines. Moreover, 
because the measured fluxes are much higher than the predictions from the accretion disk, 
self-irradiation is not likely to make a significant contribution to our measured fluxes at any 
wavelength. 
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4. SPECTRAL ENERGY DISTRIBUTION (SED) 

Using the measurements and the X-ray models in Section 2, we construct the SED of 
the ULX using the IR, optical, UV and X-ray photometric measurements. The final SED is 
plotted with the continuous black line in Figure 6. All photometric measurements described 
in the previous section were extracted from an aperture of ~50 pc centered on the ULX 
(Fig. 4). The photometric measurements presented here are likely contaminated by the 
surrounding star-forming regions west of the ULX (see the HST image in Fig. [3] and the OM 
image of Fig. 4). As a result, in constructing the SED we assume these are upper bounds to 
the ULX emission. 

X-ray to UV Band: The UV band is responsible for producing the IR emission lines, 
however it is also unfortunately the most ambiguous to fit. We found that both X-ray 
models that have power- law components extending to lower energies (PLMCD and BPL), 
over-predict the measured fiuxes if extrapolated into UV and optical. Therefore we use the 
Modified PLMCD model to extrapolate the X-ray spectrum to lower energies. Here, the disk 
emission can be approximated as a simple power- law with photon index 2/3 (or energy index 



-1/3, see lMitsuda. et al.lll984l ). since this extrapolated disk component is a superposition of 
numerous black bodies representing the integrated emission from the disk at the various 
temperatures and radii. At optical energies, this model predicts a drop in flux due to the 
finite size of the accretion disk (see Fig. 5). 

Optical Band: The disk model was further extrapolated into the optical, where the 
companion star dominates the SED. KWZ found that stars within the spectral type range 
B3Ib to 04V were consistent with optical data. Based also on the UV fluxes presented here 
we chose a B2Ib supergiant as the optical counterpart. The UV fluxes seem to be too low 
for an 05V star. The rn odel for the B2Ib star is part of the ATLAS9 Model Atmospheres 



(jCastelli fc Kurucal2004l ). with a temperature of 18500 K and a luminosity of 5x10^^ erg s~^. 
Figure 6 shows that the supergiant emits very few photons that can impact the ionization. 
On the other hand, an 05V star provides a signiflcant flux above 13.6 eV, but is equally 
deflcient above the He H edge. 

Infrared and Radio Band: At longer wavelengths, from IR to radio, we used a broken 
power-law to approximate the SED. These photometric data are viewed as upper bounds 
since the ULX flux is likely contaminated by the surrounding star-forming region. 
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5. DISCUSSION 

5.1. The SED of the ULX in Holmberg II 

In Section 4 and Figure 6 we constructed the full SED for the Holmberg II ULX using 
multi-wavelength data and previously published data. The bolometric luminosity for this 
SED is 1.34x10^° erg s~^, when the far-UV spectrum is extrapolated from the X-rays, using 
the Modified PLMCD model fit. Table 2 illustrates the degeneracy resulting from the three 
methods for X-ray spectral fitting. All three show a fairly good fit to the X-ray data, making 
it very difficult to distinguish between the different physical interpretations that each model 
represents. 

The optical and UV ffuxes plotted in Figure 6 can, in principle, come from three different 
sources: the stellar companion, the accretion disk, and from a jet. However, the overall SED 
shape and in particular the fact that the optical ffux is much lower than the UV ffuxes, 
disfavor the jet scenario as discussed below. On the other hand, the emission from the disk 
in the Modiffed PLMCD model is too faint at optical and UV wavelengths if extrapolated 
from X-rays and requires an additional source. Previously published optical data and our 
UV measurements suggest that a ~B2Ib star is the source of these ffuxes as the companion 
to the accreting BH. We explore this further in Paper 2. 

The fact that the PLMCD model overpredicts the UV and optical ffuxes if extrapolated 
to these wavelengths is expecte d, because the power-la w cannot extend too much below the 



peak of the accretion disk (e.g. iDone fc Kubotal 120061 ). On the other hand, it is intriguing 



that the simple broken power-law model (BPL) provides such a good fft to the high-quality 
X-ray data. More over, the loca t ion o f the break (~1 keV) in this case is in agreement 



with the jet model. iKaaret et al.l (120031 ) fft a very similar model to another famous ULX in 
NGC 5408. It has been shown that jets can contribute signiffcantly to the X-ray emission, 
both by optically thin synchrotron emission and Compton sc attering of the photons from 



the disk (or the companion) by relativistic electrons in the jet ( iGeorganopoulos et al.ll2002l ). 
However, even if the power-law is less steep compared to the PLMCD model, it is still 
overpredicting the UV and optical ffuxes and is therefore not the preferred scenario. 

Much better constraints are provided by photoionization modeling of the Spitzer spectra. 
The [O IV] line is only produced by high energy photons (> 54 eV), and is clearly associated 
with the ULX (Fig. [3]). For instance, the "classical" PLMCD model predicts a power-law 
component with photon index of 2.4 dominating into the UV. This model has the strongest 
emission of the three models in the far UV and soft X-rays, and should produce the strongest 
[O IV] emission. In this case, photoionization modeling can be used to predict the [O IV] 
ffux predicted from such a model, and this predicted ffux can be compared with our observed 
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values. In a companion paper (Paper 2) we present the results of such a photoionization 
analysis in an attempt to distinguish between the three X-ray models presented here. 



5.2. Characterization of the Optical/Mid-IR Line Ratios 



The ionization potential of [O IV] is just above the He II edge so both likely result from 
similar regions in the nebula and a quantitative comparison between previously published 
optical observations and the mid-IR observations can be made. KWZ found that the He II 
recombination line is consistent with photoionization by the ULX. The flux of our [O IV] 
line is almost three times higher than the He II flux measured by KWZ (2.7 x 10^^ erg s^^). 

We compared the mid-IR line ratios of the ULX with those in AGN. NGC 4395 is one 
of the few known galaxies that is likely powered by an IMBH (IFilippenko fc Hdl2003l ). It 
has a bolometric luminosity very close to our ULX, and high-excita tion emission lines were 
detected. The narrow-line He II measured by iKraemer et al. I dlgggl) f or NGC 4395 relative 
to H/3 is 0.12, a value typical for Seyfertl galaxies (jPietrich et al.ll2005l ). For the Holmberg II 
ULX, He 11/B.p ^ 0.2 (KWZ); slightly higher than for N GC 4395. The [O IV] measured from 
NGC 4395 is almost 16 times larger than the He II flux (IKraemer et al.lll999l ). However, the 
gas density is close to the critical density for this line (10^ cm~^) — much higher than the 
density of the gas surrounding the Holmberg II ULX (KWZ estimate 10 cm~^, based on the 
line surface brightness predicted with CLOUDY). The [O IV] is expected to increase as it 
approaches the critical density, while the He II is not sensitive to the gas density. An even 
more important differ ence is the abundance of O, which is ~ 5 times higher for NGC 4395 
(IKraemer et al.l 119991 ). The dependence of the [O IV] on the abundance is expected to be 
roughly linear. In spite of these two differences, we consider the [O IV] (and He II) detected 
from the Holmberg II ULX to be consistent with photoionization by a high energy source 
similar to the AGN in NGC 4395. 

The infrared [O IV] line (together with [Ne V]) is a signature of high ionization usually 



associated with actively accreting black holes (e.g. iLutz et al.lll998l ). This line is often used to 



AGN power ( 


Genze 


et al. 


1998: 


Sturm et al. 


2002; 


Satvapa 


et al. 


2004 


Smith et al. 


2004 


Satvapal et al. 


2007 




Melendez et al. 


2008 


)• 


Genzel et al. 


( 


1998) and 


Sturm et al. 


( 


2OO2I) 



have used the ratio of [O IV] /[Ne II] to distinguish between starbursts and AGN dominated 
ga laxies in the mid-IR . In Figure [7^ we compare the Holmberg II ULX with data published 
by lGenzel et al. Jl998h . The comparison clearly indicate high ionization for the ULX, similar 
to AGN. 
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Other infrared line diagnostics have been proposed by iDale et al. J200d . boosh to dis- 
entangle AGN type ionizing sources from star formation, using the ratios [S III] 33.48 nm / 



[Si II] 34.82 /im and [Ne III] / [Ne II]. Figure [7)3 shows the four regions defined by lDale et al. 



(120061 ) in their Figure 5, and their published data for different sources in the SINGS sample. 
The Holmberg II ULX is located well within region I, together with most AGNs and LINERs. 

Both diagnostic diagrams suggest that the [O IV] emission in Holmberg II is dominated 
by the ULX and that the surrounding star forming regions provide little contamination to 
this line. While this comparison provides good evidence that the line ratios for the ULX 
are very similar to AGN, we note that AGN are different than ULXs, both in their intrinsic 
emission and their environment. For example, AGN emit copiously in the UV and optical, 
and are surrounded by gas at high densities. 



Dasyra et al.l (120081 ) found a good correlation between the luminosity of the [O IV] line 
and the mass of the central BH in standard AGN. We estimate a BH mass of 3.5x10^ Mq if 
this correlation is extrapolated to the Holmberg II ULX. We can also estimate the BH mass 
using the inner disk temperature from the Modified PLMCD model (0.38 keV). We obtain a 
BH mass of 994 Mq. However, the flux of the MCD disk component does not dominate the 
X-ray emission (is only 49% of the total flux). This scenario would imply that the ULX is not 
in a typical high state, and t herefore the inner dis k temperature cannot be used reliably to 
estimate the mass of the BH (jStobbart et al.ll2006l ). Finally, using the bolometric luminosity 
calculated from our SED (1.34 x 10^° erg s~^), we estimate that the lower limit to the mass 
of the central BH, within the Eddington regime, is 85 Mq. 



6. CONCLUSIONS 

The detection of the [O IV] line from the Holmberg II ULX indicates strong UV and soft 
X-ray emission. This line is commonly used to quantify the AGN power in Seyfert galaxies. 
The ratios of this line to lower-ionization lines are comparable to those in AGNs and suggest 
that the [OIV] emission in Holmberg II results from the ULX and not from the surrounding 
star forming regions. The mid-IR line ratios indicate that the intrinsic luminosity of the 
ULX must be very high. 

The SED constructed based on the X-ray spectral modeling, UV and infrared photom- 
etry, and previously published radio and optical data suggest that the ULX companion is 
likely similar to a B2Ib supergiant. The accretion disk is not expected to contribute very 
much in the optical. The overall shape of the SED seems to favor the accretion disk emission 
as the dominant component in the soft X-ray and far UV, and is therefore likely responsible 
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for the photoionization. In the accompanying paper (Paper 2) we perform detailed pho- 
toionization modehng to provide further evidence that the detected [O IV] line is consistent 
with the ionizing emission from the ULX, and to better constrain the three X-ray models 
presented here. 
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Table 1. Measured infrared lines 



Standard Aperture Small Aperture 



Line 


Flux 


S/N Ratio 


L 


Flux 


S/N Ratio 


L 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


m 


[Ne II] 12.81 (im 


<6.08 




<0.68 


<1.98 




<0.22 


[Ne III] 15.56 fim 


12.15±3.4 


4.60 


1.35±0.39 


6.46±1.37 


7.6 


0.61±0.15 


[S III] 18.71 fim 


2.64±1.14 


4.64 


0.29±0.13 


3.02±1.16 


4.7 


0.34±0.13 


[Ne V] 24.32 


<3.88 




<0.43 


<1.16 




<0.13 


[O IV] 25.89 Mm 


7.01±1.60 


10.16 


0.78±0.2 


3.63±0.91 


44.4 


0.40±0.1 


[S III] 33.48 /im 


<4.39 




<0.49 


<1.53 




<0.17 


[Si II] 34.82 nm 


21.91±6.48 


9.90 


2.44±0.61 


6.87±1.47 


7.2 


0.65±0.16 



Note. — We show line fluxes for the standard (slightly smaller than the SH map) and the 
small aperture (4 LH pixels) as defined in Section 2. Fluxes are in 10~'^'^ W cm~^, luminosities 
(L) in 10^^ erg s""*^. If the measurements errors are smaller than the absolute calibration 
accuracy of 25%, the latter were used. For nondetections we show 3a upper limits. 



Table 2. X-ray model fits 



Model 


Nh 


kTi„/ri 


r/r2 


Ax^/dof 


L 


MOD flux 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 




(10=° cm-3) 


(keV) 






(10"° erg B-1) 




Modified PLMCD 


3.26±0.56 


0.38±0.02 


2.64±0.02 


1.02/135, 0.974/713 


1.11 


0.49 


PLMCD 


10.4±0.70 


0.27±0.a2 


2.42±0.06 


1.1386/849 


2.28 


0.12 


BPL 


5.49±0.78 


1.73±0.10 


2.57±0.02 


1.0316/849 


1.41 


0.60 



Note. — (1): X-ray model, as described in Section 2.2. (2): Intrinsic hydrogen column density. The Galactic 
column density from KWZ (3.42x10^*^ cm ^) was added separately. (3): Inner disk temperature for the MCD 
component. For the broken power-law model (BPL), this is the first (low energies) photon index parameter. 
The break is at l.U kcV. (4): Photon index for the power-law component. For the broken power-law model this 
is the second (high energies) photon index parameter. (5): Reduced values for the fit and the number of 
degrees of freedom. For the Modified PLMCD model, we quote two values for each separately fitted component 
(see text for details). (6): Unabsorbed (intrinsic) luminosities between 0.1 and 10 keV. (7): MCD component 
unabsorbed flux as fraction of the total flux for the PLMCD model. For the other two models, this simply the 
ratio of the flux between 0.1 and 1.0 keV to the flux in the whole range (0.1 — 10 keV) 
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Table 3. Multi- wavelength photometry measurements 



Instrument 


Program (observation) 


Wavelength 


Date 


Flux Density 


Aa 


Corrected Fux 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


Photometry performed in this paper 


Spitzer IRAC 


SINGS 159 


3.6 ^[m 


2004 Get 10 


0.15±0.03 m,Iy 






Spitzer IRAC 


SINGS 159 


8.0 Mm 


2004 Oct 10 


0.26±0.05 m,Iy 






Spitzer MIPS 


SINGS 159 


24.0 nm 


2004 Oet 14 


0.69±0.14 m.Iy 






XMM-Newton OM 


0200470101 


UVWl 2910A 


2004 Apr 15 


33.1±0.1 /J,7y 


0.60 


57.5 ^iJy 


GALEX NUV 


NGS 


1750-3013A 


2003 Dec 5 


26.3±0.2 pjy 


0.95 


63.1 /iJy 


GALEX FUV 


NGS 


1350-2013A 


2003 Dec 5 


33.7±0.3 (jjy 


0.84 


73.0 /iJy 


Previously published photometry 


VLA 


AT 159 


1.4 GHz 


1994 Apr 3 


1.174±0.085 mjy 






VLA 


AR165, AR227, AR258 


4.86 GHz 


1990 Jul 21 - 1991 Dec 6 


0.667±0.207 mJy 






HST ACS 


GO 9684 


P550M 


2002 Nov 24 


6.56±0.03 MJy 


0.32 


9.51 M-Ty 



Note. — All the photometric measurements performed in this paper used an aperture of 50 pc diameter (3.3 ). These values 
represent the integrated fluxes over this region, including the nebular emission and are con t amina ted by other nearby sources, 
therefore we consider them as upper bounds. The published data is taken from [Miller et al ] i2005H (radio VLA) and fro m KVZ 
(optical HST). The radio emission is extended, and is comparable to our extraction region. The optical measurement is for the 
star-like source which is believed to be the companion of the accreting BH. This measurement is likely less contaminated by the 
nebular emission or other sources. The extinctions in column 6 are both Galactic and local from KWZ. For the infrared an radio 
data the extinctions are negligible. 



(a) (b) 




Fig. 1. — Full IRS spectra obtained from matched apertures with CUBISM from the two 
apertures shown in Fig. [3] : a) The standard aperture, just slightly shorter than the SH 
aperture, b) The small aperture, this is the smallest extraction aperture possible for matched 
spectra, 4 LH pixels. By comparing the spectra it is evident that the [O IV] is clearly 
concentrated in the small aperture. 
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Fig. 2. — Details of the detected infrared lines in the spectrum of Fig. [1^. The vertical lines 
show the theoretical center for the lines at rest (dotted) and red-shifted (sohd). We used the 
recessional velocity adopted by KWZ (157 km s^^). 



[S III] 

^ [SI II] 

SI 



Fig. 3. — HST ACS false-color image of the Foot nebula in Holmberg II (left), and the 
Spitzer spectral maps of the infrared detected lines (right). The colors in the HST image 
correspond to the following filters: 814W (red), 658N (green. Ha) and 550M (blue). The 
spectral maps are: [Ne III] 15.56 /im, [S III] 18.71 //m, [O IV] 25.89 /xm, [S III] 33.48 fim and 
[Si II] 34.82 /jm. The units on the colorbar are W m~^ sr^^. The apertures used to extract 
IRS spectra (standard and the small aperture, see Section 2.1), are overplotted on all images. 
In the HST image we also show the whole LH map for comparison with the LH maps shown 
at right. The ULX optical counterpart identified by KWZ has the appearance of a bright 
star in the HST image, close to the center of the small aperture, and its position is shown 
with a green X on the spectral maps. It is located at the center of the Heel of the nebula - the 
round-shaped strong Ha emission, at the center of the small aperture. An OB association 
of stars is seen west from the ULX (see also the OM image in Fig. 4). They appear to have 
blown a large bubble in the gas, creating the "Sole" and "Toes" of the Foot nebula. The 
low-ionization lines are probably contaminated by these stars. Some contamination is also 
produced by the bright foreground star, seen on the HST image just inside the standard 
aperture, close to its northern edge. 
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Fig. 4.— Spztzer IRAC (8 /xm filter, left) and XMM-Newton OM (UVWl filter, right) images 
of the Foot nebula. The aperture used for the ULX photometry is shown as a small circle. 
At the Holmberg II distance, it has a diameter of 50 pc. With dotted lines we show the two 
apertures used to extract the IRS spectra shown in Fig. [H The larger aperture (continuous 
line) corresponds to the whole LH map. Notice that the bright star which contaminates the 
infrared emission in Fig. [3] is very faint in the UV image. The strongest UV emission in the 
nebula is located west of the ULX and it originates from the OB association better seen in 

Fig. m . 
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Photon Energy (eV) 



Fig. 5. — Examples of irradiated disks spectra. The higher flux spectrum is calculated for a 
1000 Mq intermediate-mass BH binary and an outer radius of 10^"^ cm. The lower spectrum 
is for a stellar-mass BH binary of 10 Mq and two disk sizes: 3x10^^ cm and 10^^ cm. The 
continuous lines show the total spectrum, dotted lines the spectrum from inner disk, not 
affected by the irradiation, and the dashed lines from the outer part. The vertical line marks 
the [O IV] ionization edge. The reprocessed photons have energies below 10 eV, and therefore 
not ionizing. All cases assume accretion at Eddington luminosity. 
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Fig. 6. — This fi gure shows the SED constructed in Section 4. The radio data is from 



Miller et ahl (120051 ) . the V-band magnitude for the optical counterpart of the ULX is quoted 
from KWZ. Our measured fiuxes from IRAC, MIPS, OM and GALEX are upper bounds, 
and therefore our model is below these data points. The two stellar spectra correspond to the 
spectral types range consistent with the colors and magnitudes in KWZ. The X-ray data is 
from the long (~100 ks) XMM-Newton observation ID 0200470101. The Modified PLMCD 
model was obtained by fitting the MCD and power-law components separately as explained 
in the text. The other two X-ray models, PLMCD and BPL are also shown extrapolated 
to lower energies, they clrealy overpredict the UV and optical fiuxes. The [O IV] edge at 
54.93 eV is shown vertical line. 
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Fig. 7. — Ratio line diagnostics comparison between Holmberg II ULX and other extra- 
galactic ionizing sources: AGN, starburts, etc. Its position in both plots indicate a source 
very similar to AGNs and LINERs. a) Diagnostic that uses the [ O IV] line nor r aalize d to the 
[Ne II] line. The AGN, starburts and ULIRG data is taken from lGenzel et al.l (11998). Many 
of the starburts and ULlRGs are actually upper limits (see Fig. 3 of lGenzel et al.lll998l ). For 
Holmberg II ULX we show the lower limit because [Ne II] was not detected, b) Diagnostic 
plot that uses the [S HI] and [Si II] l ir ies to distinguish between different sources. We plot 
all the data presented by iDale et al.l (120061 ) (their Fig. 5), which includes: Seyferts, LIN- 
ERs, starburt galaxies, H II regions and extranuclear star-forming regions. For Holmberg II 
ULX we show the lower and upper limits because [Ne II] and [S HI], respectively, were not 
detected. 



